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ABSTRACT 

Two recent papers (Ghez et al. 2008, Gillessen et al. 2009) have estimated the mass of and 
the distance to the massive black hole in the center of the Milky Way using stellar orbits. The 
two astrometric data sets are independent and yielded consistent results, even though the measured 
positions do not match when simply overplotting the two sets. In this letter we show that the two sets 
can be brought to excellent agreement with each other when allowing for a small offset in the definition 
of the reference frame of the two data sets. The required offsets in the coordinates and velocities of 
the origin of the reference frames are consistent with the uncertainties given in Ghez et al. (2008). 
The so combined data set allows for a moderate improvement of the statistical errors of mass of and 
distance to Sgr A*, but the overall accuracies of these numbers are dominated by systematic errors 
and the long-term calibration of the reference frame. We obtain i?o = 8.28 ± 0.15| s t a t ± 0.29| sys kpc 
and Mmbh = 4.30 ± 0.20| s t a t ± 0.30| sys x 10 6 M Q as best estimates from a multi-star fit. 
Subject headings: blackhole physics — astromctry — Galaxy: center — infrared: stars 



1. INTRODUCTION 



The motions of stars in the immediate vicinity of 
Sgr A* have been tracked since 1992 at the NTT /VLT 
and since 1995 at the Keck telescope (|Eckart &: Genzell 



119961 [Ghez et al.l[Tg 98). With the detection of accelera- 
tions (Ghez et al. 2000) and the determination of the first 



orbit (jSchodel et al.ll2002T ) these measurements provided 
firm evidence for the existence of a massive black hole 
(MBH) at the center of the Milky Way, coincident with 
the radio-source Sgr A*. The stars are used as test par- 
ticles for the gravitational potential in which they move. 
Particularly important is the bright (mjc ~ 14) star S2 
(SO-2 in the Keck nomenclature) orbiting the MBH in 
15.9 years on an ellipse with an apparent diameter of 
k, 0.2". Together with radial velocity data, the astromet- 
ric data allow for a geometric det ermination of Rn, the 
distance to the Galact ic Center (jSalim fc Gouldl 119991 : 
Eisenhauer et al. 2003). The number of orbits has in- 



creased to more than 20 since then ( 


Ghez et alJ 120051 


lEisenhauer et al. 2005: Gillessen et al.l 


2009) and in par- 



ticular the full S2 orbit has recently been used for im- 
proved determinations of th e mass and of i?o (jGhez et al.l 
\2(M [Glilessen et al.l[2009T) . 

These two astrometric data sets are independent since 
they were obtained at different telescopes and were an- 
alyzed by different teams using different tools. On the 
other hand, the radial velocity data attached to the as- 
trometric data sets largely overlap, mainly because this 
technically is straight-forward. The radial velocities re- 
fer to the local standard of rest and hence the inclu- 
sion of other data into a given data set needs no spe- 
cial care. The situation is different for astrometric data, 
for which only an approximate realiza tion of an absolut e 
reference frame currently is possible (|Reid et al.ll2007f ). 
This means that the exact definition of coordinates is a 
matter of the respective data analysis, and simply merg- 
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ing two lists of astrometric positions will fail (see fig- 
ure HJ. 

The two groups (jGhez et all 120081 : I Gillessen et al.l 
l200l come to very similar concl usions for the mass o f 
Sgr A* and R . From VLT data iGillessen et all (|2009D 
derived (when using the S2 orbit only) 



i? = 8.31±0.33| st atkpc 
Mmbh = 4.29 ±0.07| St at 

±O.34| flo M 



(1) 



where the statistical error on the mass is for a fixed dis- 
tance and the second error term is due to the statistical 
error on i?o- The latter includes already the coordinate 
system uncertainty. In addition, a systematic error of 
Ai?o « ±5 (5 kpc is present, owing to the uncertainties 
of a) how much the 2002 data (close to periastron) can 
be trusted and b) the assumption that the effective po- 
tential is Keplerian. The systematic error of i?n also 
influences Mmbh since the two quantities are correlated. 
iGhez etldl (|2008l ) obtained from the Keck data set (ne- 
glecting 2002 data, and fixing the radial velocity v z of 
the MBH to 0) 



R = 8.4 ± 0.4 kpc 
Mmbh = 4.5 ± 0.4 M Q 



(2) 



The errors are of statistical nature and the mass error in- 
cludes the uncertainty induced by the error on R . Prob- 
ably, also a systematic error should be taken into account 
here. Leaving v z free yielded somewhat larger statistical 
errors, e.g. AR = 0.6 kpc. 

The aim of this letter is to investigate the combined 
data set. We show that not only the fit results are in ex- 
cellent agreement, but also that the combined data can 
be described very well by a single Keplerian orbit. The 
such combined data set also yields an improvement in 
the statistical fit errors for the derived parameters. This 
can be understood qualitatively by inspecting the data: 
Before 2002, the Keck measurements appear superior to 
our NTT/VLT data set. At that time, our group was us- 
ing the NTT with an aperture of 3.6 m, which is notably 
smaller than 10 m for the Keck telescopes. After 2002, 
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the VLT data has comparable errors to the Keck data for 
individual data points but a much denser time sampling. 
This is owed to the location on the southern hemisphere 
of the VLT and the generous allocation of telescope time 
to the Galactic Center projects. Hence, by suitably com- 
bining the two data sets, one should obtain a data set 
which combines the respective advantages. 3 

2. DATA BASE 

The data used here is shown in figures Q] and [2] We 
resign a detailed discussion of the data, since it is well 
presented in the original works and would be beyond the 
scope of this letter. 
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Fig. 1.— Astrometric data for the star S2. Blue: NTT/ VLT 
measurements. Red: Keck measurements. The lines show the 
Keplerian orbit fits for the respective data set, which do not yield 
closed ellipses in this figure due to the motion of the point mass 
with respect to the chosen coordinate systems. The small lines 
close to the origin indicate the position of the center of mass as a 
function of time. 

The astrometric data from the NTT /VLT we are us - 
ing are identical to the data from [Gillessen ct al. (2009) 
up to August 2008. Since then, one more epoch has 
been observed in 2008 and 11 more epochs in 2009. We 
tr eated these data in pre cisely the same way as described 
in iGillessen et al.l (|2009f ) , where the somewhat cumber- 
some derivation of astrometric positions from the imag- 
ing data is described in detail. As done previously, we 
assigned lower weights to the VLT data from 2002 to ac- 
count for the fact that the star might have been confused 

3 We use data collected between 1992 and 2009 at the European 
Southern Observatory, both on Paranal and La Silla, Chile; in 
particular from the Large Programs 179.B-0261 and 183.B-0100. 



in that year. In total, the NTT/ VLT data set contains 70 
epochs from March 1992 to August 2009. Furthermore, 
we added one more epoch (March 2009) to the measure- 
ments of reference stars by which the coordinate system 
is tied to the I CRF 

iGhez et"aTI (|2008l ) presented in their table 3 the mea- 
sured positions for all Keck epochs 4 . We use the given 
26 epochs, ranging between 1995 and 2008. 

For the radial velociti e s, we use the combined da ta 
base from lGillessen et al.l (|2009[ ) and lGhez et"aL1 (|2008t ) 5 . 
Finally, we add one more VLT epoch from 2009, which 
brings the total number of radial velocity points to 36. 
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Fig. 2. — Radial Velocity data for the star S2. Red: Keck mea- 
surements. Blue: VLT measurements. The black line shows al- 
ready the result of a combined orbit fit. 



3. ORBITAL FITS 

Figure Q] shows that the two astrometric data sets can- 
not be put together in a simple fashion. To account 
for the difference we make the following, simple assump- 
tion: The two data sets only differ in their definition of 
the coordinate system, namely the position of the ori- 
gin and the zero velocity 6 . This yields four parameters 
(Ax, Ay, Av x , Av y ) by which the difference can be de- 
scribed. Given that both data sets contain many more 
data points, it is feasible to solve for these four param- 
eters during the fitting procedure. These four parame- 
ters come hence in addition to the six orbital elements 
and the parameters describing the gravitational potential 
(Mmbh, Rq, position, zero velocity, v z ). 

We have im plemented the additio nal four parameters 
in the code of Gillessen ct al. (2009). For (most of) the 
orbital fits, we used the same assum ptions as for the 
preferred fit in IGillessen et alj (|2009D . namely we im- 
posed priors on the NTT/ VLT coordinate system def- 
inition which reflect our best estimate of how well an 
absolute coordinate system was established: 

V(x) = 0± 1.0 mas V{v x ) = ± O.lmas/yr 
V{y) = 0± 2.5 mas V{v v ) = ± 0.1 mas/yr (3) 
V{vJ) =0±5km/s . 

4 The definition of their X-coordinate (increasing values to the 
east), however, does not match the signs of the numbers in the 
table, which have to be reverted (S2 was West of Sgr A* in the 
1990s and was East of Sgr A* from 2003 to 2007). 

5 In table 4 of the latter paper two radial velocity measurements 
from the year 2002 are missing. They are shown in the figures 
and clearly have b een used in the fit. The values can be found in 
IGhez et aj] <2003ft though. 

6 A difference in rotation is very unlikely, since the images can 
be oriented very well by the means of stars across the field of view. 
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Using that, the combined data set can be described very 
well by a Keplerian orbit (figure [3]). For the 211 degrees 
of freedom (96 astrometric epochs which count twice and 
36 radial velocity points minus 17 parameters) we ob- 
tained x 2 = 271 for the fit shown in figure [3l correspond- 
ing to a reduced \ 2 of 1.28. 
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Fig. 3. — Result of the combined orbit for the star S2. Blue: 
NTT/VLT measurements. Red: Keck measurements. The black 
line shows the Keplerian fit (row 1 in table |2j- 

The fit result for the four additional parameters is 

Ax = —3.7 ± 0.6 mas 

Ay = -4.1 ±0.6 mas 
Av x = -0.68 ± 0.11 mas/yr 
Av y = 0.26 ± 0.10 mas/yr . (4) 

We conclude that these parameters are well-determined 
from the fit, which validates our approach for combining 
the two data sets. The positional difference (Ax, Ay) 
is somewhat larger than the actual mismatch between 
the two coordinate systems. This happens simply since 
(Ax, Ay) refer to the epoch May 2005, the zero point 
in time of the NTT/VLT coordinate system, while the 
Keck data uses July 2004 for that. Nevertheless it is 
worth checking how well the numbers in equation [3] com- 
pare to the expected accuracies of the coordinate sys- 
tems. The uncertainty of the NTT/VLT coordinate sys- 
tem (equation[51) is smaller than the values in equation^) 
IGhez etldl (|2003l ) conclude that their coordinate system 
is accurate to (Ax, Ay) = ±(5.7, 5.7) mas in position 
and to (Av x , Av y ) = ±(0.6, 0.9) mas/yr. This is consis- 
tent with equation [4j The larger uncertainty of the Keck 
coordinate system is due to the shorter timeline used for 



TABLE 1 

Best fitting orbital elements for S2 from a combined 
orbit fit. The errors are rescaled such that the reduced 
X 2 is 1. The corresponding gravitational potential is 

GIVEN IN THE FIRST ROW OF TABLE [2] 

Parameter Value 

~~ a 0.1246" ± 0.0019" 
e 0.8831 ± 0.0034 

i 134.87° ± 0.78° 

SI 226.53° ± 0.72° 

ui 64.98° ±0.81° 

t P 2002.3293 ± 0.0066 



determining the motions of the reference stars (« 2 years 
compared to w 6 years for the NTT/VLT data). 

In tabic [T] we give the orbital elements for the com- 
bined fit. The corresponding parameters describing the 
gravitational potential are presented in the first row of 
table [21 We continued by making different choices for 
the priors, for the 2002 VLT data and for the model of 
the gravitational potential. These fits are also listed in 
tabled the respective choices are indicat ed in each row. 

Com paring the results with table 4 of iGillesscn et all 
(2009) one notices that generally the statistical fit errors 
for Mmbh and Rq are 20% to 30% smaller in the com- 
bined orbit fit. Also the errors on the position and 2D ve- 
locity of the MBH have decreased moderately. Generally, 
the fitted position of the MBH has moved by « 2.5 mas 
East and « 1.5 mas North. This is consistent with the 
accuracy by which the radio source Sgr A* can be iden- 
tified in the infrared images. 

The comparison also shows that for the combined data 
set the various fits differ less from each other than for 
the NTT/VLT-only data set. For example, including or 
not including the 2002 VLT data (comparing the fi rst 
and second row) changed Rq in lGillessen et al.1 (|2009f ) by 
0.95 kpc, while in the combined data set the same change 
is only 0.62 kpc. Also the sixth row (a fit excluding the 
2002 VLT data and neglecting priors) sh ows that the re- 
sults a re much closer in the combined set. Gillcsse n et al.l 
(2009) reported Rq — 6.63±0.91 kpc, while the combined 
data yield Ro = 7.73 ± 0.57 kpc. The range of values 
for Rq is in the com bined data set 1 kpc only, while for 
IGillessen etHI (<2009D it exceeded 2 kpc. 

Next, we turned to post-Newtonian orbit models. A 
relativistic orbit fit is shown in row 9 of table [21 yield- 
ing a slight increase in Rq compa red to the Kepleria n 
fit. This bias was already noted bv lZucker et a l. (2006). 
Ass uming the relativis tic f3 2 effects for the radial veloc- 
ity (|Zucker et a l. 2006) we then checked if the combined 
data allow to constrain the relativistic periastron shift. 
For this purpose we introduced a parameter /, which is 
for a Keplerian model and 1 for the correct relativis- 
tic model. The fit yielded a very large error A/ w 10 
and we conclude that also the combined data set is not 
yet sensitive to the relativistic precession induced by the 
Schwarzschild metric. 

We also repeated the tests of Gillcss en et al.l (|2009[ ) 
for an extended mass component, using three different 
power law slopes for the density profile: p(r) ~ r~ a , 
a = 1.4, 7/4, 2.1. The normalization of the density pro- 
file was a parameter, which we determined by the fit. We 
express the results for this parameter in table [2] by giving 



4 



Gillessen et al. 



TABLE 2 

The cen tral potential from fit ting the S2 data. The first six rows use the same assumptions as the first six rows in 

TABLE 4 OF lGlLLESSEN ET AL.I (120091) AND THUS ARE DIRECTLY COMPARABLE. THE LAST FOUR ROWS SHOW POST-NEWTONIAN FITS. ALL 
ERRORS ARE RESCALED SUCH THAT THE REDUCED X 2 IS 1. THE SECOND COLUMN IS A FLAG INDICATING WHETHER THE 2002 VLT DATA 
FROM S2 WAS USED. THE THIRD COLUMN INFORMS WHICH OF THE PRIORS FROM EQUATION [3] WERE USED. COLUMN FOUR IS A FLAG 
WHETHER A RELATIVISTIC OR A KEPLERIAN MODEL WAS USED. THE FIFTH COLUMN INDICATES THE RADIAL DENSITY PROFILE OF THE 
ASSUMED EXTENDED MASS. THE LAST COLUMN IS THE FITTED MASS BETWEEN PERI- AND APOCENTER OF THE S2 ORBIT DIVIDED BY MyiBH- 





VLT 


priors 


rcl. 


ext. 


Ro 


AfMBH 


a 


<5 


v a 


vs 


Vz 


n 




2002 






mass 


(kpc) 


(10 6 Af o ) 


(mas) 


(mas) 


(^as/yr) 


(fias/yr) 


(km/s) 


% 


1 


yes 


2D, v z 


no 




8.34 ± 0.27 


4.40 ± 0.27 


2.93 ± 0.55 


0.41 ± 0.85 


-61 ± 84 


89 ± 91 


2.5 ± 6.3 




2 


no 


2D, v z 


no 




7.72 ± 0.33 


3.90 ± 0.29 


3.21 ± 0.58 


-1.53 ± 1.13 


-109 ± 88 


67 ± 94 


0.4 ± 6.5 




3 


yes 


v z 


no 




8.39 ± 0.28 


4.44 ± 0.29 


3.18 ± 0.62 


0.14 ± 0.88 


-27 ± 115 


186 ± 132 


2.2 ± 6.2 




4 


no 


Vz 


no 




7.67 ± 0.37 


3.87 ± 0.31 


3.72 ± 0.66 


-2.05 ± 1.21 


-134 ± 123 


106 ± 138 


0.24 ± 6.4 




5 


yes 


none 


no 




8.77 ± 0.41 


4.99 ± 0.53 


3.13 ± 0.62 


0.41 ± 0.88 


-34 ± 115 


173 ± 131 


33 ± 24 




6 


no 


none 


no 




7.73 ± 0.57 


3.93 ± 0.58 


3.69 ± 0.67 


-1.95 ± 1.37 


-131 ± 123 


107 ± 138 


3.7 ± 27.0 




7 


yes 


2D 


no 




8.78 ± 0.41 


5.03 ± 0.55 


2.87 ± 0.55 


0.68 ± 0.85 


-63 ± 84 


81 ± 91 


37 ± 24 




8 


no 


2D 


no 




7.85 ± 0.52 


4.04 ± 0.55 


3.19 ± 0.58 


-1.34 ± 1.23 


-109 ± 88 


66 ± 94 


7.9 ± 26.7 




9 


yes 


2D, v z 


yes 




8.46 ± 0.28 


4.65 ± 0.30 


2.58 ± 0.55 


-0.10 ± 0.84 


-37 ± 84 


97 ± 91 


1.5 ± 6.3 




10 


yes 


2D, v z 


yes 


r" 1 ' 4 


8.55 ± 0.29 


4.70 ± 0.31 


3.03 ± 0.66 


0.25 ± 0.87 


-130 ± 114 


70 ± 94 


1.4 ± 6.3 


2.0 ± 1.6 


11 


yes 


2D, v z 


yes 


r -7/4 


8.51 ± 0.28 


4.66 ± 0.30 


3.07 ± 0.66 


0.14 ± 0.86 


-140 ± 114 


67 ± 93 


1.1 ± 6.3 


1.8 ± 1.4 


12 


yes 


2D, v z 


yes 


r- 21 


8.49 ± 0.28 


4.64 ± 0.30 


3.09 ± 0.66 


0.00 ± 0.85 


-146 ± 115 


64 ± 93 


0.6 ± 6.3 


1.6 ± 1.3 



?7, the mass between peri- and apocenter of the S2 orbit 
divided by the fitted mass of the MBH. The results in 
rows 9-12 are not dire ctly comparable to the numbers in 
IGillessen et al.l (|2009l ) who imposed an additional prior 
on i?o (obtained from the other S-stars) for the fits with 
an extended mass component. The results from the com- 
bined data yield nevertheless very similar upper limits: 
At most very few percent of the central mass can make 
up an extended mass component inside the S2 orbit. 

Overall, the fits presented here show that the chosen 
scheme to combine the two independent data sets is valid 
and that consequently (moderately) smaller statistical 
errors are obtained. Still, the approach is not yet suffi- 
cient to detect post-Newtonian effects in the S2 orbit. 

4. DISCUSSION 

Since more than 10 years two groups have assembled 
independently astrometric data sets for the stars in the 
central arcsecond of the Milky Way. The data sets are 
truly independent. They were obtained at different tele- 
scopes with different instruments. The analysis chains 
used different tools (deconvolution in the NTT/ VLT 
case, PSF-fitting for the Keck data), and the definition 
of astrometric coordinates is imple mented in d i fferent 
way s. While it was reass uring, that iGhez et al.l (|2008l ) 
and IGillessen et al.l (|2009[ ) concluded very similarly for 
]1/mbh and Rq , it was not clear that the agreement actu- 
ally also holds for the underlying data. We were able to 
show that indeed the most simple assumption - namely 
that the two data sets only differ in the coordinate sys- 
tem definition - is sufficient to perfectly map the two data 
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sets on top of each other. 

We presented for the first time a combined orbit fit, 
which however only moderately improves the accuracy by 
which Mmbh and R§ can be derived from the S2 data: 

i?o = 8.34± 0.27| sta t ± 0.5| sys A:pc 
Mmbh = 4.40 ± 0.27| stat ± 0.5| sys x 10 6 M . (5) 

For the sake of completeness, we cite also the updated 
numbers for a fit using not just the combined S2 data, 
but in addition SI, S8 , S12, S13 and S14, the same stars 
IGillessen et~aTI (|2009D had used: 

Ro = 8.28 ± 0.15| stat ± 0.29| sys kpc 
Mmbh - 4.30 ± 0.20| stat ± 0.30| sys x 10 6 M & . (6) 

Besides the expected improvement in Ai?o| s tat from 
0.17 kpc to 0.15 kpc, also the systematic errors have de- 
creased mildly compared to the previous work due to the 
smaller influence of the S2 2002 data. 

We conclude that the combination does not help much 
in overcoming current limitations. The true value of the 
two data sets actually is their independence, allowing for 
cross checks and lending credibility to the results. 

Further substantial improvements in measuring the 
gravitational potential from Sgr A* by means of stellar 
orbits probably will come from improved astrometry on 
existing data, from longer time lines with the existing in- 
struments and finally from advances in instrumentation. 
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